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ABSTRACT 

The issue of giant planet formation by core instability (CI) far from the central star is rather contro- 
versial because the growth of massive solid core necessary for triggering the CI can take longer than 
the lifetime of the protoplanetary disk. In this work we assess the range of separations at which the 
CI may operate by (1) allowing for arbitrary (physically meaningful) rate of planetesimal accretion 
by the core and (2) properly taking into account the dependence of the critical mass for the CI on the 
planetesimal accretion luminosity. This self-consistent approach distinguishes our work from similar 
studies in which only a specific planetesimal accretion regime was explored and/or the critical core 
mass was fixed at some arbitrary level. We demonstrate that the largest separation at which the CI 
can occur within 3 Myr corresponds to the surface density of solids in the disk > 0.1 g cm~^ and 
is 40 — 50 AU in the minimum mass Solar nebula. This limiting separation is achieved when the 
planetesimal accretion proceeds at the fastest possible rate, even though the high associated accretion 
luminosity increases the critical core mass delaying the onset of the CI. Our constraints are indepen- 
dent of the mass of the central star and vary only weakly with the core density and its atmospheric 
opacity. We also discuss various factors which can strengthen or weaken our limits on the operation 
of the CI. 

Subject headings: planets and satellites: formation — protoplanetary disks — stars: planetary systems 



1. INTRODUCTION. 

Recent discoveries of distant planetary companions 
(separations of at least tens of AU) around several nearby 
stars by direct imaging (Marois et al. 2008; Thalmann et 
al. 2009) have stimulated investigation of their formation 
mechanisms. Two contending theories of the planet for- 
mation — the core instability (Perri & Cameron 1974; 
Harris 1978; Mizuno 1980; hereafter CI) and the gravi- 
tational instability (Cameron 1978; Boss 1998) — differ 
quite dramatically in their ability to form planets at var- 
ious distances from the star. 

It is generally thought that CI should be capable of 
forming giant planets close to the star, at separations 
< 10 AU. Both theoretical modelling (Mizuno 1980; 
Stevenson 1982; Pollack et al. 1996) and the existence of 
two gas giants at semi-major axes of 5.2 AU and 9.5 AU 
in our own Solar System (coupled with some knowledge 
about the evolution of the planetary orbital architecture 
in the early Solar System) attest to this statement. At 
the same time it seems unlikely that giant planets can be 
formed by the gravitational instability in the inner parts 
of protoplanetary disks because of the long local cool- 
ing time of gas (Matzner & Levin 2005; Rafikov 2005, 
2007). On the contrary, giant planet formation by the 
direct gravitational instability appears feasible at large 
distances from the star (Boley 2009; Rafikov 2009; Clarke 
2009; but see Boss (2006) for an opposite view) where the 
disk cooling time is short enough to permit fragmenta- 
tion of the gravitationally unstable disk (Gammie 2001). 
Whether the CI can operate beyond several tens of AU 
from the star is not so clear. 

The well known problem faced by the CI far from the 
star is that the buildup of massive refractory core (neces- 
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sary for triggering the vigorous gas accretion) is thought 
to take very long time, longer than the 1 — 10 Myr life- 
time of the gaseous component of protoplanetary neb- 
ula. Studies illustrating this problem are usually based 
on two key assumptions (Dodson- Robinson et al. 2009): 
first, that the CI is triggered whenever a Merit ~ 10 M^ 
core is built by planetesimal accretion, and, second, that 
the planetesimals accrete onto the growing core at rather 
modest rates (Ida & Lin 2004; we will explain later what 
do we mean by that). When these conditions are imposed 
it is generally found that the CI does not commence at 
a > 20 — 30 AU prior to the nebular gas removal. 

In this note we show that aforementioned assumptions 
are too restrictive and neither of them need to be adopted 
in determining the feasibility of the CI. First, planetes- 
imal accretion by the growing core can (at least poten- 
tially) be much faster than what has been previously as- 
sumed. Second, the critical core mass itself strongly de- 
pends on the planetesimal accretion rate. Based on these 
observations we formulate in this work a novel constraint 
on the operation of the CI in protoplanetary disks. 



2. PLANETESIMAL ACCRETION. 

Onset of the CI is intimately related to the accretion of 
solid material by the protoplanetary core, as discussed in 
more detail in [JS] For this reason we start by reviewing 
the process of planetesimal accretion by growing cores. 

It is well known that the rate at which a core ac- 
cretes planetesimals is a strong function of planetesi- 
mal random velocities (Dones & Tremaine 1993; here- 
after DT93). An important parameter for determining 
the dynamical state of planetesimal population is the ra- 
tio p = Rc/Rh of the core radius Rc to the Hill radius 
Rh = a{Mc/Mi,Y/'^ , where a and Mc are the semi-major 
axis and the mass of the core, and M^, is the mass of the 
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star. At large separations from the star 

is much less than unity, where pc is the core density and 
we define am = a/(10 AU), M^^ = M^,/Mq, p^.i = p/1 
g cm~'^. 

DT93 have demonstrated that whenever p ^ 1 there 
are four possible regimes of planetesimal accretion. 

• When the planetesimal random velocity disper- 
sion^ a > nRHV^^'"^ = y/GMjc^ is the an- 
gular frequency at the core's location) - a regime 
that we call very high dispersion - accretion is very 
slow because gravitational focussing is inefficient 
and the collision cross-section is given by the geo- 
metric cross-section of the core. 

• When QRh ^ c ^ Q,RhP~^^^ {high dispersion in 
the notation of DT93) the relative planetesimals 
velocity with respect to the core Vrei ~ cr. Grav- 
itational focussing is important and planetesimal 
accretion rate M grows as a goes down. 

• When flRnp^^'^ ^ flRn {intermediate disper- 
sion) one finds that Vrei is no longer determined by 
the random motions of planetesimals but is rather 
set by the shear in the differentially rotating disk, 
so that Vrei ^ i^Rn- Gravitational focussing satu- 
rates at a constant value whenever a < flRn and 
M keeps increasing as cr decreases simply because 
the thickness of the disk goes down increasing local 
volume density of accreting bodies. 

• Finally, when a < VIRhP^^^ {very low dispersion) 
disk becomes so thin that a core can accrete the 
whole vertical column of material doomed for col- 
lision with it. Given that focussing is at its max- 
imum this situation corresponds to the most effi- 
cient accretion regime. 

The rate of accretion in the very low dispersion case is 
(DT93) 

M = Mraax ~ (6.47 ± Qm)np^ '""^ sRr . (2) 
« 5.1 Gi/^M^/^E.p-i/e^ (3) 

where Es is the surface density of solids. In this study 
we parametrize as 

= lSiaro"gcm-', (4) 

where Si = S^/ (1 g cm^^) and a is a constant (minimum 
mass Solar nebula (MMSN) corresponds to a w 3/2 and 
El « 1). Note that Mmax is independent of M^,. 

For comparison, it is often assumed that protoplane- 
tary cores predominantly accrete planetesimals with the 
dispersion of random velocities a ~ ^IRh, which corre- 
sponds to the transition between the so called shear- and 
dispersion-dominated dynamical regimes (e.g. Dodson- 
Robinson 2009). This is what we called a "modest" ac- 
cretion rate before. At this transition DT93 find 

Mtr « (6 - 10)pJ:s^R%, (5) 

^ For simplicity here we do not differentiate between the hori- 
zontal and vertical velocity dispersion. 



and one immediately sees that Mtr ^ Mmax since p 1. 

Previous studies of the CI have been limited to con- 
sidering planetesimal accretion only in the high- velocity 
regime (cr > flRa) or at the boundary with the interme- 
diate velocity regime (i.e. for a w HRh)- While large 
(tens to hundreds km in size) planetesimals indeed likely 
accrete in these dynamical regimes at rather slow rates 
it is quite possible that most of the mass growth of the 
core is not due to accretion of these large bodies. Accord- 
ing to Rafikov (2004a) and Goldreich et al. (2004), as the 
core becomes massive enough it dynamically excites large 
planetesimals in its vicinity since damping agents such as 
gas drag are not effective for massive bodies. This leads 
to fragmentation (rather than growth) of planetesimals 
when they collide with each other. Resulting fragmenta- 
tion cascade converts significant fraction of the solid mass 
into small mass debris. These fragments finally reach 
small enough sizes that they are dynamically "cooled" 
by gas drag to low velocities, which makes it possible 
for them to be accreted at very high rate compatible 
with Mmax (Rafikov 2004a, b) in the very low dispersion 
regime. Recent coagulation calculations done by Kenyon 
& Bromley (2009) confirm this general picture, lending 
support to the possibility of maximally efficient accretion 
in the very low dispersion regime. 

In this work we are interested in obtaining a robust 
limit on the operation of core instability, which means 
that we need to consider all possible modes of plan- 
etesimal accretion and actually determine the accretion 
regime facilitating the onset of the CI the most (instead 
of selecting it ad hoc). It is not obvious a priori which 
accretion regime is best for achieving the CI within the 
limited amount of time: on one hand high M allows to 
build the core of a given mass faster, but on the other 
hand high M also implies larger critical mass (see Q, 
which takes longer to build. 

We will discover in SjHthat fastest route to the CI lies 
through the most efficient accretion at the rate Mmax 
and for this reason we parametrize planetesimal accretion 
rate as 

/ M \ 

M = xMmax ~ 6A7xp'^^^sna^ f , (6) 

«2 X lO-^Me yr-ixanTSiP.^i/'M.'/^o, 

{Mc^io = Afc/10 M^) where the parameter x accounts 
for the deviation of M from Mmax- In general x may be 
a function of Mc and a, and it is expected that x ^ 1 
(although in some situations X ^ 1 may possible, see 
fj6]). But we will show in ^ that slow accretion (i.e. 
X ^ 1) results in smaller distance from the star at which 
planets could form by the CI than in the case x ^ 1- 
It should also be emphasized that our concentration on 
the highest M accretion regime is not only for the sake 
of the argument - as mentioned before this regime may 
actually occur quite naturally in the course of the core 
buildup. 

Accretion luminosity corresponding to the M given by 
equation ^ is 

L = « 8.2xS,G3/2 Af4/3pi/6 (7) 

Rc 

«8xlO-*^L0 x«m"Sip'/fA<fo. 
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According to equation ([6]) accretion at the rate M^ax 
results in the core radius growing hnearly with time^. 
Time needed for the core mass to reach a predetermined 
value Mc by accretion at the rate x^max is [see equation 
@] 

1/6,^1/3 1/6 ,.1/3 

r(M,) « 0.6-^-^ « 0.3 Myr a^^—^—. (8) 

It follows from this formula that a core can grow to 10 
in 3 Myr by accretion at the rate M^ax even at 
~ 50 AU in the MMSN. Previous results which did not 
account for the factors leading to the fast accretion (plan- 
etesimal fragmentation and gas drag on debris) predict 
much smaller distance (< 20 AU) at which such core 
would be able to form on that time scale. 

Having said all that, in the next section we show that 
the formation of the 10 core is not a prerequisite for 
a CI and that the Merit is a function of accretion rate. 

3. PROTOPLANETARY ATMOSPHERE. 

Numerical calculations (Mizuno 1980; Pollack et al. 
1996; Ikoma et al. 2000) and analytical theory (Stevenson 
1982) show that CI commences when the core mass Merit 
is so large that the mass of its atmosphere is comparable 
to the core mass itself, 

Matm{Mcr^t) = VMcrtt, (9) 

i.e. when the gaseous component of the protoplanet be- 
comes self-gravitating. Here 77 is a parameter of order 
unity; its exact value may depend on the accretion his- 
tory of the core (Ikoma et al. 2000). Condition ^ should 
be regarded as an equation for Merit which can be eas- 
ily solved as long as the dependence of Matm on Mc 
and other parameters of the problem is specified. Thus, 
calculation of Merit involves understanding properties of 
planetary atmospheres and calculation of Matm in par- 
ticular. 

It has been known since the works of Mizuno (1980) 
and Stevenson (1982) that the critical core mass is a 
rather weak function of the density po and temperature 
Tq of the surrounding nebula (see Rafikov [2006] and 
^5 5.31 for the discussion of the conditions under which 
this behavior is expected). This has led to the wide ac- 
ceptance of the CI idea since the interiors of the Solar 
System giant planets are believed to harbor ~ 10 M^ 
cores despite their different separations from the Sun. 
Unfortunately, this observation of the Merit invariance 
has shadowed the fact first noticed by Stevenson (1982) 
that Merit is a strong function of the planetesimal accre- 
tion luminosity L, or, equivalently, planetesimal accre- 
tion rate M. This fact has been subsequently confirmed 
both numerically (Ikoma et al. 2000; Hori & Ikoma 2010) 
and analytically (Rafikov 2006). 

Stevenson's analytical (1982) results strictly apply only 
to cores possessing radiative atmospheres with constant 
opacity. Rafikov (2006; hereafter R06) has studied more 
general types of atmospheres and has shown that they 
segregate into two classes depending on whether L is 

One can easily see that accretion at the rate Mtr (see Eq. [5]) 
1 /3 

also leads to oc t, although the coefBcient of this relation 

is considerably smaller than in the case of M = Mmax and has 
different scaling with a. 



higher or lower than some critical luminosity Lcr- When- 
ever L > Lcr the intense energy release near the core sur- 
face makes protoplanetary atmosphere and the nebular 
gas in the Hill sphere of the core convectively unstable. 
Such atmosphere have high entropy and rather low mass 
relative to the mass of the core. In the opposite case, 
when L < Lcr, protoplanetary atmosphere is separated 
from the nebular gas by a roughly isothermal shell of gas 
in which energy is transported radiatively and gas en- 
tropy decreases from the nebular value to a much smaller 
value characteristic for the atmosphere. Density in this 
shell increases roughly exponentially towards the planet 
which makes Matm niuch higher (for a given core mass) 
than in the high-luminosity case. 

It was demonstrated in R06 that far from the star pro- 
toplanetary atmospheres are virtually always character- 
ized by L ^ Lcr, even if planetesimal accretion proceeds 
at the maximally efhcient rate Mmax- For that reason 
we will consider only the low luminosity atmospheres in 
this work. The total atmospheric mass Matm in the low 
luminosity case was computed in R06 by making simpli- 
fying assumptions about the conditions in the deep layers 
of the planetary atmosphere: either a polytropic model 
with constant polytropic index (mimicking the fully con- 
vective interior) or a fully radiative atmosphere with a 
simple power-law parametrization of the opacity depen- 
dence on gas pressure and temperature. Under these as- 
sumptions a significant fraction of the atmospheric mass 
resides in the outer^ layer, near the inner boundary of 
the roughly isothermal external radiative zone, in agree- 
ment with calculations by Mizuno (1980). It can then be 
shown that (R06) 

where is the Boltzmann constant and ^ is a weakly 
varying factor, which can be computed exactly for a given 
density distribution inside the atmosphere. This formula 
shows how Matm depends on important parameters like 
Mc and kq and on accretion history, to which L is sensi- 
tive. 

However it is not obvious that the simplifying assump- 
tions about the atmospheric properties used in deriving 
formula pop are justified given the complexity of the 
physical effects encountered deep in the atmosphere: self- 
gravity gradually becomes important as Mc approaches 
Merit, grain sublimation changes opacity in a non-trivial 
fashion, the equation of state may be varying with depth 
because of molecular dissociation, an atmosphere may 
have both radiative and convective regions at the same 
time. These complications likely do not affect the qual- 
itative conclusions reached in R06 but should be impor- 
tant for quantitative comparisons. 

On the other hand, numerical calculations of the CI 
which include the aforementioned physical effects (Ikoma 
et al. 2000) are typically limited to exploring the depen- 
dence of Merit on only a limited set of input parameters 
such as kq and M. The latter is usually taken to be con- 
stant through the calculation and this assumption signif- 
icantly limits the direct applicability of these numerical 

^ This is different from the case of an atmosphere with constant 
opacity considered by Stevenson (1982). 
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results to realistic situations (since typically M increases 
as Mc grows), including our present study. 

To be able to apply the existing numerical results for 
the cases when the core accretion history is non-trivial 
{M ^ const) while at the same time keeping the flex- 
ibility of the theory outlined in R06 we have resorted 
to the following approach. First, we still use formula 
(fTO|) to calculate Matm but now we do not assume co- 
efficient C to be constant as various physical effects in 
the deep atmosphere introduce additional variations of 
the density profile not captured by the analytical theory 
of R06. Second, using equation © we calculate Merit 
under the assumption M = const used in Ikoma et al. 
(2000). With this information in hand we can calibrate 
the dependence of C on various physical parameters of the 
problem by comparing our results with those of Ikoma et 
al. (2000) . Finally, after C has been calibrated for a par- 
ticular accretion history M =const, equation (fT0| can be 

applied for a more general situation, e.g. for the Af (Mc) 
dependence used in this work. 

This approach should work well as long as the state of 
the atmosphere is fully determined by the current value 
of M, and is independent of the full accretion history. 
For this to be the case the thermal timescale of the at- 
mosphere must be short compared to the planetesimal 
accretion timescale, and this has been verified in R06. 
This gives us confidence that the outlined method should 
be robust and justifies its application to the problem at 
hand. 

Since C is predominantly affected by the changes of 
the thermodynamical state of material in the atmosphere 
(grain sublimation, molecular dissociation, variations of 
the equation of state depend on gas temperature and 
density) one should expect ( to depend most strongly on 
Mc and L — variables that determine T and p deep in 
the envelope. Ikoma et al. (2000) have shown that in 
their case Merit scales with M and kq roughly as power 
laws. This motivates us to look for the dependence of ( 
on Mc and L also in the power law form. The details of 
the calibration procedure are presented in Appendix |X] 
where we show that 



Co « 4 X 10 



(11) 



-32 „.-i5 



82 M, 



-<5 



provides a reasonably good fit to the numerical results 
of Ikoma et al. (2000). Note that C, has been found to 
depend on L only weakly, and can be considered a func- 
tion of Mc only. In deriving this result for C we have 
assumed that the coefficient rj appearing in the relation 
^ is independent of the planetesimal accretion history 
of the core and is thus the same in Ikoma et al. work as 
well as in our case. 

4. CRITICAL CORE MASS. 

Results obtained in the previous section allow us to 
compute Merit for M given by equation ([6]) : substituting 
Matm in the form (ITUl) with ( given by equation (fTTT) into 
the instability condition ^ and using equation (O for 
core's accretion luminosity we find 

- „ X 4 „ 1/6 13/(5+35) 



Mc. 



(12) 



^40 Mfl 



xSiKo.ip;l^f 
[A 



0.35 



-0.35a 



where kq.i = ko/(0.1 cm^ g~^) and p2 = /^/(2m//) (rn^ 
is the atomic hydrogen mass). 

Equation (IT^ shows among other things that Merit is 
indeed a function of the planetesimal accretion regime: 
according to equation ^ we may think of the free pa- 
rameter X as a direct measure of planetesimal M, and 
Merit scales roughly as ^■'^^ , Just for illustration, if (as 
has been done in Dodson-Robinson et al. 2009) planetes- 
imal random velocities were kept at the level a ^ flRn 
(rather than being essentially zero as assumed in deriv- 
ing equation ([6])) the equation ([6|) could still be used but 
with X « P^^'^ < 1 (DT93; R06). At a = 40 AU where 
p « 10~^ one would then find Merit smaller by a factor 
of 5 compared to what equation (|12|) predicts for x = I 
at the same location. 

With equation we also show for our chosen M 
behavior that Merit varies with the distance from the 
central object: for x==const the critical mass scales as 
« fl-o-s^ i-jjg MMSN (a = 3/2). This makes Mertt at 
10 AU twice as large as Merit at 40 AU, everything else 
being equal. Thus, taking Merit to be independent of the 
regime of planetesimal accretion and distance from the 
central object is generally not justified. 

5. LIMIT ON CORE INSTABILITY. 

For the CI to happen before the nebular gas dispersal 
the growth time of the core with mass Merit must be 
less than the lifetime of the protoplanetary nebula T„eb, 
i.e. T{Merit) < Tneb- Plugging in our result ([T2|) into 
equation ([8]) this constraint can be rephrased in terms of 
the lower limit on the surface density of solids: 



{2+S)/2 , ^4 



-| 1/(4+35) 



^'J neb ^ 



(13) 
(14) 



iO.l g sm 



3 Myr J 



-1.13 „0. 21^0. 13 
Pel '^O.l 



This inequality represents the main result of this work — 
a robust lower limit on the planetesimal surface density 
at which CI is capable of producing giant planets within 
a protoplanetary nebula lifetime. Note that Ejim is a 
sensitive function of the nebula lifetime Tneb, while it 
depends rather weakly on the bulk density of the core 
Pc, atmospheric dust opacity kq, and the mean molecular 
weight of the atmospheric gas p,. 

Now we can also determine the limiting core mass Mum 
defined as the critical core mass for x^s = ^Hm, i-e. at 
the very extreme of the region where CI can still occur 
within Tneb- It is found by substituting into ([T^ : 



Miim — 



4^ 



P- 



4 1/3 
Pc KO 



3/(4+3.5) 



(15) 



-0.4 ^0.13^0.4 
PcA '^O.l 



3 Myr / 



By construction Mum is independent of the surface den- 
sity profile in the nebula. This mass is to be compared 
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with the isolation mass Mi^o (the core mass at which it 
has accreted all sohd mass within its feeding zone) — an 
annulus centered on core's orbit and having a full width 
equal to S,Rh' 



M- — 



,2^ 



3/2 



1/2 



(16) 



'.2 {^/5f'^'/'a%'Ms 



(an MMSN-like density profile was used in making nu- 
merical estimate). At 44 AU one finds Miso « 6 M^, 
which is smaller than Mum- However, a modest radial 
displacement of the core due to some type of migration 
can easily expose it to additional fresh material allowing 
Mc to reach Mum (Alibert et al. 2005). Alternatively, 
increasing Eg (boosting up Ei) by a factor of 2 makes 
Miso ~ Mum at 44 AU. 

Quite interestingly, the value of Mum is not too 
far from 10 commonly accepted as the core mass 
throughout the whole protoplanetary disk. This coinci- 
dence is accidental since our estimate of Mum was de- 
rived self-consistently rather than postulated in an ad 
hoc fashion. 

5.1. Limiting distance for core instability. 

Given a constraint ()14|) and having a particular model 
of the radial distribution of one can determine the 
maximal radial extent aum of the region in protoplan- 
etary disk in which the CI can produce giant planets 
within the nebula life time r„e&- Since is expected to 
be a decreasing function of a this would only be possible 
for a < aiim- Using our power- law parametrization (j4]| 
of E^ wc find that in a MMSN-likc disk with a = 3/2 



.ffjf'^ « 44 AU 



(XS] 



^2/3 0.35 
) P2 



„0. 14^0. 09 



^0.1 



3 Myr 



0.75 



(17) 



Accretion at M — M^^ax corresponds to x = 1 in this 
formula. The limiting distance found in equation (jl7p is 
similar to previous estimates (e.g. Dodson-Robinson et 
al. 2009) obtained for less vigorous planetesimal accre- 
tion (i.e. for smaller M) and fixed Merit = 10 but 
this is just a coincidence. 

Equation (|17l) shows that as Es (or Ei) increases, the 
extent of CTcapable part of the protoplanetary disk also 
grows. But Eg cannot be increased without limit — at 
some point the gaseous component of the protoplane- 
tary nebula would become self-gravitating, and, depend- 
ing on its cooling time (Gammie 2000), would either 
fragment or evolve quasi-viscously while maintaining the 
marginally gravitationally unstable state (Rafikov 2009; 
Clarke 2009). Thus the most extreme value ofaum can be 
obtained by taking Es(a) — fdgT,Q^i{a) in equation ([T4|) . 
where fdg is the dust to gas ratio, which we take to be 
10^^, and EQ=i(a) is the gas surface density at which the 
disk is marginally gravitationally unstable, which hap- 
pens when the Toomre Q = Qcsfi/{TrG'E,Q^i) is of order 
unity (Safronov 1960; Toomre 1964). Assuming that the 
disk is heated by a central star with luminosity 1 Lq at 
normal incidence (i.e. neglecting complications related 
to flaring geometry, Chiang & Goldreich [1997]) we find 
that 



— 7/4 

fdg^Q=i{a) « 20a^o g cm 



(18) 
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< 
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0.1 



10 
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Fig. 1. — Plot of the limiting semi-major axis au^ beyond 
which core instability cannot produce planets within 3 Myr as a 
function of Si — surface density of solids at 10 AU in a MMSN- 
like protoplanetary disk with Ss oc a~^/^ (see Eq. [4]). Solid curve 
is for planetesimal accretion proceeding at maximum possible rate 
Mmax, dotted curve is for M = Mtr (accretion of planetesimals 
with (T ~ flR[f, see Eq. [5]). Dashed lines correspond to in 
a marginally gravitationally unstable disk, both for M = Mmax 
(upper line) and M = Mtr (lower line). See text for details. 



Using this density profile in equation (1141) we find 

0.64 



a 



Q=l 
lim 



200 AU 



,^0.57, ,0.3 
A M2 
„0. 12^,0. 07 
Pel '^0.1 



3 Myr J 



(19) 



In practice this limit will hardly apply to real systems be- 
cause self-gravitating disks with Q ~ 1 would not persist 
for the lifetime of the nebula, as we already mentioned. 

5.2. Sensitivity to planetesimal accretion efficiency. 

Through our calculations we have retained in all formu- 
lae the parameter x defined in equation ([5]), which char- 
acterizes the efficiency of planetesimal accretion. This 
allows us to see how the variation of planetesimal M 
with respect to Mmax affects the planet formation by 
the CI. From equation we see that smaller x (corre- 
sponding to less efficient accretion) results in higher Eg 
and smaller aum, see equations ([T7)) . (fTO|) . 

As we emphasized in fj3l this result is not trivial since 
less efficient accretion implies lower planetesimal lumi- 
nosity, bigger Matm for the same Mc, and lower Merit i 
making CI easier to get going. However, it turns out that 
the competing effect of being able to grow the massive 
core faster at higher Af is more important for the CI to 
commence within the limited amount of time. For this 
reason the largest extent of the region in which the CI 
can happen in time T„e6 is reached when the core is able 
to accrete at the highest possible rate, namely at Mrnax 
(we discuss whether it is potentially possible to exceed 
Mmax in ^6 6.2p . This makes our numerical estimates in 
equations (ITT)) . (fT9)) very robust. 
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Just for illustration let us also consider a situation^ in 
which the core accretes planetesimals at the rate Mtr 
given by equation ([5]), and determine aum in this case. 
Comparing expressions ([2]) and ([5]) one can easily see 
that we can do this by simply setting x = P^^^ {pi'^) 
is defined by equation [I]) either in equation (jl4p or in 
equation (I17p and determining aum from the resulting 
expression: 

where we assumed MMSN-like disk properties to provide 
direct comparison with affJl^'^^ given by equation (|T7)) . 
If, as we did in deriving equation (jl9|) , we assume that the 
disk is marginally gravitationally unstable we would find 
^?irn^'*^ « 20 AU for our standard choice of parameters. 
In Figure [1] we display different expressions for aum as 
functions of Si — value of at 10 AU — obtained for 
Ko.i = 1, Pel = 1, /i2 = 1, Tneb = 3 Myr under a variety 
of assumptions regarding and x- 

Clearly, less efficient accretion results in significantly 
more compact region of the disk where the CI can oc- 
cur. In fact, equation (EUl) implies that the formation of 
Saturn by the CI at 10 AU in our Solar System would 
only be possible if the surface density in the proto-Solar 
nebula was at least 4 times higher than in the MMSN or 
if the nebula dissipation timescale was at least 10 Myr, 
in agreement with existing studies which assume plan- 
etesimal accretion at M ^ Mtr (Ida & Lin 2004). 

5.3. Sensitivity to opacity variations. 

Rather interestingly, our calculations find very weak 
dependence of both E;j„j and a/^^ on kq. Although the 
reduction of kq does help to reduce Merit and Mum, 
the values of Ej^m and aum remain virtually unaffected. 
Given the large uncertainly in the value of kq this prop- 
erty further strengthens our estimates of Ejim and aum- 

Recently Hori & Ikoma (2010) calculated Merit as a 
function of M for protoplanets with dust-free (possible 
if dust grains sediment from the outer layers of the pro- 
toplanetary atmosphere) and metal-free (i.e. containing 
only H and He) atmospheres and found values of Merit 
lower by up to an order of magnitude (as low as 2 Mq in 

the metal- free case for M ~ 10"^ Mq yr~^) than in the 
dusty case. It would certainly be interesting to repeat 
calculations done in §ij4l[5]and Appendix 1X1 for the case 
of dust-free atmosphere to see how this extreme reduc- 
tion of opacity would extend the radial range available 
for CI. 

In practice, we cannot do such calculation at present 
since it is not possible to calibrate Matm against the 
results of Hori & Ikoma (2010) as we did in equations 
(fTO)) -([TT |) and Appendix [X\ This is because in the ab- 
sence of dust /t is a function not only of gas tempera- 
ture but also of gas density. As demonstrated in R06 
in this case Merit is no longer independent of the ambi- 
ent temperature Tq and density po of the nebula as the 
original analyses of Mizuno (1980) and Stevenson (1982) 
suggest. Instead one finds that Matm oc {po/Ti)''/'^^+''\ 
where q = d\nK/d\np (R06). Obviously, Merit then 
also depend on po and Tq, and the knowledge of this de- 
pendence is very important for obtaining 'Sum and aum 



in the dust-free case. Unfortunately, we do not possess 
this knowledge from first principles as opacity calcula- 
tions are rather complicated, and in any case we cannot 
currently calibrate Merit as functions of Me, Tq and po 
against numerical results because calculations of Hori & 
Ikoma (2010) were done for a single value of the plane- 
tary semi- major axis (meaning fixed values of Tq and pq), 

while M was varied. The scaling of Matm and Merit with 
Tq and pq and its implication for the possibility of the CI 
thus remain worthwhile issues for future investigation. 

We can still get a qualitative idea of how aum changes 
in the dust-free case by setting opacity in equation (fT7| 
at the very low level consistent with pure gas opacity, e.g. 
K = W)-^ cm2 g-i. We then finda,f^=f^^ « 80 AU com- 
pared to 44 AU that equation (jl7| predicts for k = 0.1 
cm^ g""'^- Thus, opacity reduction due to sedimentation 
and coagulation of dust grains in the protoplanetary at- 
mosphere may help in extending the range of distances 
in which the CI is possible. 

6. DISCUSSION. 

Despite the robustness of our arguments it is not incon- 
ceivable that some additional factors can weaken them 
and make giant planet formation by the CI possible even 
beyond the limits represented by equations and (IT7)) . 
Alternatively, it is quite possible that some of the as- 
sumptions used in deriving these results are too extreme 
and one can get even better constraints by focussing on 
less dramatic assumptions. Below we review factors that 
can work one way or another. 

6.1. Extending CI to larger radii. 

One possible way to facilitate CI and increase aum is 
to consider possibility of planetesimal accretion at rates 
exceeding Mmax- This is very difficult (since there are 
many factors that tend to reduce M compared to Mmax, 
see §6 6.2|) but may be possible if e.g. one takes into ac- 
count the increase of planetesimal capture cross-section 
by the core caused by its extended, dense atmosphere. 
This effect has been previously investigated by Inaba & 
Ikoma (2003) who demonstrated that an increase of M 
by a factor of ~ 10 compared to the value computed with- 
out atmosphere is possible. According to equation (|17p 
such an enhancement of M (incorporated by increasing 
X) would boost a^Jf^^ by a factor of 4. 

Our present calculations assume that the core is accret- 
ing planetesimals continuously until the protoplanetary 
nebula dissipates — this is important at large a since 
massive core requires long time to be built. But one 
may wonder if building smaller core in shorter time and 
then cutting off subsequent planetesimal accretion (and 
energy release at the core surface, which supports atmo- 
sphere against going unstable) completely may still lead 
to the CI and potentially extend it to larger semi-major 
axes. Such accretion scenario has been adopted by e.g. 
Pollack et al. (1996). The problem in this case is that 
even if Af = it still takes long time for the atmosphere 
around the core to grow to the mass comparable to Me. 
INEO show that this process occurs on thermal timescale 
of the atmosphere and typically takes millions of years. 

Similar problem is also encountered in a scenario where 
the core grows rapidly by planetesimal accretion in the 
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inner regions of protoplanetary disk and then gets scat- 
tered out to large radii by some massive perturber. One 
might expect that after the orbit of the scattered core 
circularizes by dynamical friction the core would gradu- 
ally accrete massive atmosphere and undergo CI at some 
point. Given that both the orbit circularization and en- 
velope accretion are likely to take long time it is not at all 
obvious whether the CI could be achieved in this scenario 
within several Myrs. 

6.2. Limiting CI to smaller radii. 

There are many factors that can potentially reduce 
aiim compared to 44 AU estimated in equation (|17p . For 
example, there are several reasons why it may not be 
possible for M to reach the maximum rate ^rnax • 

First, the growing core may clear out a gap in planetes- 
imal disk around its orbit, thus significantly reducing M 
(Tanaka & Ida 1997; Rafikov 2001; 2003a). In our previ- 
ous calculations we implicitly assumed this not to hap- 
pen e.g. because of the core migration through the disk, 
which allows fresh planetesimal material to be constantly 
supplied for core accretion (Alibert et al. 2005). 

Second, as we mentioned in SQl a known pathway to 
M « Mmax is via the growth of the core to the size at 
which it starts dominating dynamical evolution of nearby 
planetesimals and triggers their efficient collisional frag- 
mentation (Rafikov 2004). However, there is a strong 
implicit assumption in this scenario — that the core can 
reach this critical size within the nebula lifetime. Rafikov 
(2003b) has shown that at a ~ 30 — 40 AU a dynami- 
cally dominant core would need to have mass of order 
10^"* g and would require on the order of 10 — 100 Myr 
to grow in the MMSN. This time scale is apparently in 
conflict with the typical dissipation times of protoplane- 
tary disks. Thus, one may need to either require a more 
massive disk at these radii or to find other pathways for 
accretion at Mmax- Formation of massive solid bodies 
by direct gravitational instability facilitated by various 
streaming instabilities (Johansen et al. 2009) may be 
quite relevant for the latter option. Regardless of this 
(arguably rather serious) issue our estimate pT|) still re- 
mains a useful upper limit on aum- 

Among other factors hindering the onset of CI and re- 
ducing aiim we should mention the possibility of high 
opacity in the protoplanetary atmosphere. It was sug- 
gested in ij5 5.31 that k may be very low because of the 
dust sedimentation and growth. However, infalling plan- 
etesimals which feed core accretion likely get partly dis- 
rupted in the atmosphere leaving behind large amount of 
refractory material. This may actually increase k com- 
pared to the value of 0.1 cm^ g~^ assumed in equation 
(fTT]). Nevertheless, given the weak sensitivity of aum to k 
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the potential increase of atmospheric opacity is unlikely 
to have huge effect on aum- 

7. SUMMARY. 

We studied the formation of giant planets by core in- 
stability at different locations in the protoplanetary disk 
with the goal of determining the range of radii where 
the CI is feasible within the several Myr lifetime of the 
protoplanetary disk. We demonstrate that this range is 
determined by two key factors: 

• The high planetesimal accretion rate is necessary to 
build the solid core as rapidly as possible at large 
separations from the star. 

• Intense energy release at the core surface caused 
by planetesimal accretion increases the critical core 
mass and delays the CI. 

The first factor turns out to be more important and the 
largest distance at which CI can happen, around 40-50 
AU, is obtained when M is at its highest possible value 
corresponding to accretion of dynamically cold planetesi- 
mals. The core mass corresponding to this case is around 
15 Mg, likely compatible with the isolation mass at this 
distance. 

Our approach is quite different from other similar stud- 
ies which often assume that (1) accretion proceeds at 
much slower rate Mtr defined by equation ([5|) correspond- 
ing to accretion of planetesimals moving with random 
velocities at the level of ilRn and/or (2) the CI com- 
mences after the core has reached a fixed mass of around 
10 irrespective of the planetesimal accretion rate or 
the location in the disk. Relaxing these two arbitrary as- 
sumptions we are able to obtain a significantly more ro- 
bust and self-consistent limit on the CI operation which 
can be represented as a lower bound on the solid surface 
density (> 0.1 g cm~^) or an upper bound on the size 
of the region where the CI can get going within several 
Myr timescale. These limits are insensitive to the mass 
of the central star and depend only weakly on the opac- 
ity in the core atmosphere. Our predictions are relevant 
for interpreting the results of current and future direct 
imaging surveys (Marois et al. 2008) designed to uncover 
and characterize the population of giant planets at large 
separations from their parent stars. 

Author is grateful to Ruth Murray-Clay for useful dis- 
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The financial support for this work is provided by the 
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APPENDIX 
CALIBRATION OF C- 

INEO investigated the dependence of Merit on kq and M assuming the latter to be constant. Their numerical 
calculations of the CI include processes such as the dissociation and ionization of hydrogen, realistic gas opacities 
with the inclusion of the dust grain contribution, and self-gravity of the atmosphere. These are the crucial ingredients 
missing in the analytical models of Stevenson (1982) and R06. Thus, it is expected that calculations of INEO should 
be more realistic than the aforementioned analytical studies. 

We calibrate the coefficient C entering the expression for Matm against the results of INEO who found that for 
kq > 10^^ cm^ g~^ (see equation A2 in INEO) 



7Mff 



M 



10-7 Me yr- 



Kq 



1 cm^ g 



and it follows from their Table 6 that the best-fit power law exponents q' and s satisfy q' 
Motivated by these scalings we first assume power law dependence in the form 



(Al) 



0.27-0.29 and s = 0.24-0.29. 



C = Cor?MfL-". (A2) 
The explicit dependence of ( on parameter 77 is motivated by the core instability condition ([9]) and the assumed 
independence of 77 on the planetesimal accretion history (i.e. the same value of 77 applies in our case of M cx Mc'" , see 



eq. [6], as in the M ~ 
equation (|9]) we find 



where 



:const case studied by INEO). Plugging this expression into ()10|) and using L — GMcM / Rc and 



q = 



3(1 + 0;) 
7 + 3(5 - 2a; ' 



(l+w)/3' 



7 + 3(5-2w' 



(A3) 



(A4) 



play the role of g' and s in INEO case. It follows from (jA4[) that qjs = 1 + while the INEO's results for these 
indices give 1 < g'/s < 1.2. This constrains cj < 0.2 resulting in a very weak dependence of C, on L. Based on this 
observation and given the approximate nature of our calibration procedure we decided to neglect the dependence of 
C on L altogether and to consider ^ in the form dTTI) . i.e. scaling with only. Repeating our calculation we find 
instead of (|A3[) the following expression for Merit- 



Merit — 



1/3- 



s = q = 



7 + 3(5 



(A5) 



Choosing 8 — 1.2 as stated in equation (|lll) we find s ^ q - 
INEO. 

To determine Co we use the fact that according to equation (jAl[) M^^ 
Me yr"-'^. This uniquely determines the value of Co in equation (fill) . 

Note that in carrying out this calibration we implicitly used the fact that Merit is independent of the ambient 
conditions in the nebula if k is independent of gas density (R06), which is true if opacity is independent of gas density. 
This is an important point since the INEO calculations have been done without varying the external conditions. More 
on this issue can be found in §5 5.31 



0.28 in good agreement with the q' and s values found in 
i 7 Me for = 1 cm^ g^^ and M = 10^'^ 



